t=0:0.01:5;
sys = tf([1],[1,3,2]);

f = exp(-2*t).*heaviside(t);
r = lsim(sys,f,t);

plot(t,r);

